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Abstract. Preliminary results from HI and ZEUS on Deep Inelastic Scattering (DIS) 
at high momentum transfer squared Q 2 are presented. Used are all available e + p data 
accumulated by the HI and ZEUS experiments between 1994 and 1997, corresponding 
to integrated luminosities of 37pb _1 and 47pb _1 , respectively. The anomalies observed 
at high Q 2 in the 1994 to 1996 data still remain, though with less significance. Since 
this high Q 2 domain represents a new frontier in DIS, the same data are used to 
search for new particles possessing direct couplings to lepton-quark pairs. Assuming 
that the slight excess of events observed in Neutral Current DIS is due to a statistical 
fluctuation, preliminary limits on the production of leptoquarks and of squarks in 
R-parity violating MSSM are presented. 
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1. Introduction 

Deep Inelastic Scattering (DIS) of leptons by nucleons at the ep collider HERA offers 
an unique possibility to probe the proton at very small distances (~ 10~ 16 cm). The 
domain of very high momentum transfer squared Q 2 is also a window where first signs 
of physics beyond the Standard Model (SM) may manifest themselves. 

Using 1994 to 1996 data, corresponding to a total integrated luminosity of about 
35pb~\ HI and ZEUS [[| experiments reported an observation of an excess of 
events in DIS at Q 2 > 15000 GeV 2 . This finding stimulated a lot of experimental and 
theoretical activity. Data accumulated during 1997 allowed to increase the statistics by 
roughly a factor 2.5. The measurements of DIS at high Q 2 based on all available e + p 
data are presented and compared to the SM expectations in section |2]. The same data 
are then used to search for new particles possessing direct couplings to lepton-quark 
pairs, namely leptoquarks and squarks in SUSY models where R-parity is violated ($ p ). 
Preliminary results on leptoquark production are shown in section |3], together with 
some remarks on particular phenomenological aspects at HERA. Analogous constraints 
on resonant ^-squarks production are presented in section |j. 

2. Deep Inelastic Scattering at high Q 2 at HERA 

2.1. Comparison of DIS data with the SM expectations 

The full available e + p data have been used by HI and ZEUS experiments to update 
the measurements of Neutral Current (NC) and Charged Current (CC) DIS presented 
in H |J. The distribution in Q 2 of the NC DIS events measured by HI |J is shown 
in Fig. 0a (all data, £ = 37.4pb -1 ) and Fig. 03 (1997 data alone, C = 22.8pb~ x ), in 
comparison to the SM expectation. Also shown in Fig. |l]b and Fig. |]d are the ratios 
of these measurements to theory. The curves above and below unity specify the ±1<t 
levels determined from the combination of statistical and systematic errors of the NC 
DIS expectation. These errors are dominated by the uncertainty on the electromagnetic 
energy scale of the calorimeter and vary between 8.5% at low Q 2 and about 30% at the 
highest Q 2 . The NC DIS expectation agrees well with the data at Q 2 < 10000 GeV 2 , 
while at larger Q 2 deviations are observed with a slight deficit of events between 
10000 GeV 2 and 15000 GeV 2 and a slight excess at Q 2 > 15000 GeV 2 . The 1997 data 
alone suggest similar deviations, although with a more marginal significance. Table [l] 
summarizes the numbers of events observed by both experiments and expected in the 
SM with Q 2 | above various thresholds Q 2 min . At Q 2 > 35000 GeV 2 , ZEUS observes two 
events while 0.29 are expected. These events were reported previously @. In the 1997 
data, no additional events were found in this Q 2 domain. 

| Q\ and Q\ a mean Q 2 reconstructed by the electron and double angle method, respectively. 
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Figure 1. (a) Q 2 dis- 
tribution of the observed 
NC DIS events in com- 
parison to the SM expec- 
tation (histogram); (b) 
ratio of the observed 
and expected number of 
events as a function of 
Q 2 ; (c,d) same as (a,b) 
but for the 1997 data 
alone. 
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Table 1. Num- 
ber of observed (N a t s ) 
and expected (N exp ) NC 
DIS events with momen- 
tum transfer squared Q 2 
above given thresholds 



o 2 

mm 


Exp. 


N obs 


1 ' exp 


Q 2 e > 15000 GeV 2 


HI 


22 


14.8 ±2.1 


Q 2 a > 15000 GeV 2 


ZEUS 


20 


17 ±2 


Q 2 e > 25000 GeV 2 


HI 


6 


1.6 ±0.3 


Q 2 a > 35000 GeV 2 


ZEUS 


2 


0.29 ±0.02 



2.2. Differential Cross Sections 

Differential cross sections in terms of Q 2 are derived by both experiments |3|, f|, 
For that purpose, the number of events in each Q 2 bin is corrected for acceptance and 
migration, using Monte Carlo simulation, and for photon radiation and electroweak 
effects, using the HERACLES || program. Fig. § shows the NC and CC cross sections 
da/dQ 2 obtained by HI and ZEUS. 

For NC DIS, the rapid fall off of da/dQ 2 by about 6 orders of magnitude is on 
the whole well described by the SM prediction. However, as already mentioned, the 
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HERA e + p DIS cross section 94 - 97 



Figure 2. Dif- 
ferential cross sections 
da/dQ 2 for NC and 
CC DIS, measured by 
HI (squares) and ZEUS 
(points), using all avail- 
able e + p data. The 
curves represent the SM 
predictions. 
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Figure 3. Left: ratio of the NC DIS cross section measured by HI to the SM 
prediction. Right: ratio of the CC DIS cross section measured by ZEUS to the SM 
prediction. On both figures, the inner bars represent the statistical error and the outer 
ones the total error. 



measurements slightly undershoot the SM expectation at Q 2 ~ 10000 GeV 2 , while 
a slight excess of events is observed in the two highest Q 2 bins. The ratio of the 
experimental (HI) and theoretical NC cross section is shown in Fig. |3] (left). Here, the 
SM prediction is obtained from a NLO QCD fit performed by HI || , which includes the 
new HI high Q 2 data as well as HI and fixed target charged lepton-nucleon data at lower 
Q 2 . The inner error bars on the data points correspond to the statistical errors, while 
the outer ones indicate the full error, the typical total systematic error being about 4%. 

For CC DIS, the Standard Model also provides a very good description of the data 
over most of the kinematic range as can be seen in Fig. 0. Only at the highest Q 2 the 
data show a tendency to lie above the SM prediction. The ratio of the measured (ZEUS) 
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to the theoretical (CTEQ4) da/dQ 2 distribution for CC DIS is shown in Fig. § (right). 
The shaded region indicates the theoretical uncertainty due to the uncertainties in the 
partons densities. As before, the outer error bars on the data points correspond to the 
total error. The main source of systematic error is the calorimeter energy scale. It leads 
to an error on the measured cross section of 5 — 10% in most of the Q 2 range, except at 
the highest Q 2 , where the error rises to about 50% H. 

The CC differential cross section da/dx has also been measured by ZEUS and 
compared to the SM prediction. One can see some tendency in the data to lie above 
the expectation. If the CTEQ4 parametrization of parton densities is modified to give a 
larger d/u ratio at high x [0], the excess observed at high x relative to CTEQ4 becomes 
less significant [§]. 



3. Leptoquark Searches 

3.1. Phenomenology 

Leptoquarks (LQ) are scalar or vector bosons which carry both lepton and quark 
quantum numbers. They appear in many extensions of the SM. At HERA, LQs may 
be produced resonantly, via the fusion of the incident lepton with a quark or antiquark 
coming from the proton. 

A reasonable phenomenological framework is provided by the BRW model |J. This 
model is based on the most general Lagrangian that is invariant under SU (3) x SU(2) x 
£7(1), respects lepton and baryon number conservation, and incorporates dimensionless 
couplings of LQs to left- and/or right-handed fermions. Under these assumptions LQs 
can be classified into 10 isospin multiplets, half of which carry fermion number F = 
and couple to e + + q, while the other half carries F = 2 and couples to e + + q. The 
F = (F = 2) class contains 2 (3) scalar and 3 (2) vector multiplets. In order to 
avoid stringent bounds from low energy experiments which would put LQs outside the 
reach of present day colliders, couplings are allowed to only one combination of fermion 
chiralities, and only within a single family. The couplings are generically called A and 
are unknown parameters of the model. In the above model, first generation LQs only 
decay into e + q and v e + q with fixed branching ratios, (3{eq) + (3{vq) = 1. We will use 
the nomenclature of [H| to label the different LQ states. 



In the following analysis, it is assumed that only one of the LQ multiplets is 
produced at a time, and that the states within a given isospin multiplet are degenerate 
in mass. At HERA, LQs can be resonantly produced in the s-channel and exchanged in 
the u-channel. This is illustrated by the diagrams in Fig. Obviously, the s- and it- 
channel processes do not interfere with each other, but each of these channels interferes 
with the SM 7 and Z exchange. The relevant matrix elements can be found in f{|||. 

§ The signs of the interference terms between standard model gauge boson contributions and fermion 
number F = leptoquark contributions given in |^] are incorrect. The correct signs can be inferred 
from and amount to multiplying the right-hand-side of Eqs (14d,e) and (17d,e) by (-1). 
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Figure 4. Illustra- 
tive diagrams for reso- 
nant LQ production in 
the s-channel (a) and 
u-channel LQ exchange 
(b). 




LQ 



(a) 





(b) 



Figure 5. Ratio of the 
cross section err for a— 
channel LQ production 
and finite LQ width T 
to the narrow-width ap- 
proximation. 
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For the allowed strength of the couplings and in the mass range interesting for 
HERA, the decay width of a LQ into a lepton and a quark is generally small : 
T s = 3/2IV = X 2 M/16n ~ 40 MeV for M = 200 GeV and A = 0.1, S and V referring 
to scalar and vector states, respectively. In the narrow- width approximation (NWA), 
the resonant production cross section is given by o~nwa = f;A 2 g(x), q(%) being the 
probability to find the relevant quark q in the proton with a fraction x = M 2 /s of the 
proton momentum. In the above, y/s ~ 300 GeV is the collision energy in the ep centre 
of mass and M is the LQ mass. 

For high mass LQs, i.e. at high x, the parton densities in the proton are falling 
very steeply. At the same time, the LQ width increases with M leading to important 
tails at smaller values of x. As a result, the convolution of the parton density with the 
corresponding Breit-Wigner distribution gives a considerable larger cross section than 
the NWA. This is quantified in Fig. [5|, for the typical range of couplings which can be 
probed at HERA with the current luminosity. One can see that for F = and A = 0.5 
the deviation from NWA becomes noticeable at M > 250 GeV, while for F = 2 the 
finite width becomes non-negligible already at M > 200 GeV. 

Furthermore, the u-channel contribution can also become large especially for F = 2 
LQs and e + p scattering where valence quarks are involved. A third effect which plays 
an important role at high M is the interference with SM 7 and Z exchange. Depending 
on the LQ type, the interference can be constructive or destructive, and it can get 
comparable to or greater than the pure s- and u-channel contributions. In Fig. || the 
relative magnitudes of the various terms in the high Q 2 domain are shown for the Sq,r 
(F = 2) and the Si/ 2) l (F = 0) scalar leptoquarks. For the So,r, the w-channel and 
the (positive) interference with SM DIS begin to dominate the s-channel cross section 
at M > 220 GeV. For the Sx/2,l, the w-channel contribution is negligible, but the 
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Figure 6. Con- 
tributions of the s-and 
u-channel processes and 
of their interference with 
SM DIS for (a) S ,r 
(F = 2) and (b) S 1/2<L 
(F = 0). A has been 
fixed to 0.5, a value 
which is typical for the 
experimental sensitivity 
in this mass range. 
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(negative) interference becomes important for M > 285 GeV. 



3.2. Analysis 

Resonant leptoquarks decaying into an electron and a quark jet lead to signatures which 
resemble NC DIS events. However, the signal could be identified as a resonance peak in 
the invariant mass distribution, with events in the peak showing an angular distribution 
specific to the spin of the leptoquark state. For example, the isotropic decay of a 
scalar LQ in its CM frame leads to a flat spectrum in y, where y = (1 + cos#*)/2, 
6* being the polar angle of the scattered lepton in this frame. This is remarkably 
different from the steeply falling 1/y 2 distribution expected at fixed x for the dominant 
t-channel 7 exchange in NC DIS events. Hence, a LQ signal should become more and 
more prominent at high y (or high Q 2 = yM 2 ). Consequently, it should be possible to 
optimize the signal to NC DIS background ratio by a mass dependent lower cut on y 
(or equivalently on cos 8*). Note that this optimal cut depends not only on the spin 
of the LQ, but also on its fermion number, since the s- and w-channels, which lead to 
different angular distributions, contribute differently for F = and F = 2. 

The NC DIS candidates at Q 2 > 2500 GeV 2 observed by HI are shown in Fig. 0a 
in the (M, y)-plane. Here, the full dataset was re-analysed including a new in-situ 
energy calibration for the electrons. This has led to slight migrations (within originally 
quoted systematic errors) of individual events in (M,y). Fig. ^b shows the projected 
mass spectra, before and after applying the optimal y cut mentioned above. In the mass 
window considered in U of total width 25 GeV around 200 GeV, and at y > 0.4, 8 events 
are observed, while 3.0 ±0.5 are expected from NC DIS. Out of these, 5 events originate 
from the 1994 — 1996 data and 3 from the 1997 data. Hence, the "mass clustering" 
around ~ 200 GeV observed previously is not confirmed by the 1997 HI data alone. 

ZEUS searched for a resonance in positron-jet mass M e j, without applying any 
constraints. The M e j does not give the optimal resolution for LQ, however, it is free 
from initial state QCD and QED radiation shifts. The mean reconstructed mass is 
determined to be within 6% of the generated value, whereas the peak is on average 4% 
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Figure 7. High Q 2 NC DIS candidates observed by the HI (a,b) and ZEUS (c,d) 
experiments, (a) Events in the (M, y)-plane. Contours of fixed Q 2 are shown, (b) Mass 
spectra of the events before (triangles, dotted lines) and after (circles, solid lines) the 
optimal y cut. (c) Events in the (cos 9*, M ej )-plane. Shaded is the removed region 
which corresponds to the interface between calorimeters. Open circles show the high- a; 
and high-y events reported previously, (d) Mass spectra of the events. 



lower than the one obtained from a LQ MC. Fig. [7]c shows the distribution of ZEUS NC 
DIS candidates in the (M, cos6 l *)-plane. The fiducial cut applied on the positron which 
excludes the interface region between the central and forward calorimeters is indicated 
as the shadowed area. The projected spectrum of the invariant mass M e j is shown in 
Fig. ^d. It is interesting to note that the outstanding high Q 2 events discussed in [[J 
tend to cluster around M ej ~ 215 GeV. At M ej > 200 GeV and cos6»* > -0.2, 4 events 
are observed which is in slight excess to the 1.2 events expected from NC DIS. On the 
other hand, for the events at M e j > 200 GeV the angular distribution of the scattered 



positron is found to be similar to that expected from the SM [13 
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Figure 8. Left : HI exclusion limits at 95% CL on the Yukawa coupling A as a 
function of the LQ mass (a) for F = 2 and (b) for F = scalar leptoquarks. Right : 
analogous limits from ZEUS (a) for scalars with F = 2, (b) vectors with F = 2, (c) 
scalars with F = 0, and (d) vectors with F = 0. 



3.3. Results 

Assuming that the slight excess of NC events observed at high Q 2 is due to a statistical 
fluctuation, HI and ZEUS have derived preliminary limits on the Yukawa coupling A as 
a function of the LQ mass, for all LQ species classified by BRW . These are shown in 
Fig. ^. The discrepancy between the HI and ZEUS bounds for high mass scalar LQs is 
due to the use of the NWA approximation by ZEUS in the calculation of the signal cross 
section, while HI takes into account the LQ width by a Breit-Wigner distribution. For 
Yukawa couplings of electromagnetic strength, A 2 /47r = cvem, the existence of F = 
LQs is excluded up to M ~ 275 GeV. 

Relaxing the BRW assumptions, the HI Collaboration also obtained bounds in 
the plane (3 vs M, (3 being the branching ratio for LQ — > eg, for fixed values of the 
Yukawa coupling. The results can be found in [12] and in these proceedings [|14]]. For 



(3 = 1, TeVatron experiments exclude first generation LQ masses below ~ 242 GeV |15[ . 
independently of any assumption on A. Nevertheless, for low /3, an important discovery 
potential remains for HERA. For example, if j3 = 0.1, the D0 experiment excludes 
masses only up to about 110 GeV, while HI rules out masses below 255 GeV if A = 0.3, 



and 210 GeV if A = 0.1 [0. The DELPHI p| and OPAL |Tj experiments at LEP 
recently performed a search for single LQ production and exclude LQ masses below 
142 GeV for A = 0.3, which because of the smaller centre of mass energy cannot be 
competitive with the sensitivity achieved at HERA. 
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Figure 9. Bounds on X^j, 

j = 1,2, against the £1/2, l mass, 
for several fixed values of An. 
The grey domains are excluded at 
95% confidence level. The full 
curve indicates the most stringent 
indirect low energy limit on A31 for 
A n = 0.3. 
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HI also carried out a search for Lepton Flavor Violating LQs looking for 
LQs decaying into /i + gj]] and r + qj through the coupling X 2 j and X 3 j respectively. 
Mass dependent bounds on X%j are shown in Fig. |9] for several fixed values of An. For 
An = A 3j = 0.3, HI excludes LQs lighter than 255 GeV, which significantly improves 
the limit of about 100 GeV obtained from third generation LQ searches carried out at 
the TeVatron |Tj|, [2T][] , unless A is much smaller. Indirect limits on A31 are also improved 
by typically one order of magnitude. 



4. Searches for R-parity violating squarks 



4-1. Phenomenology 

The most general superpotential consistent with the gauge symmetry and supersymme- 
try of the Minimal Supersymmetric Standard Model (MSSM) contains Yukawa terms of 
superfields which violate lepton (L)- and baryon (5)-number. These interactions can 
give rise to rapid proton decay. In order to avoid this problem, one introduces additional 
discrete symmetries which forbid some or all of the dangerous Yukawa couplings. A rad- 
ical solution is achieved by requiring R-parity conservation where R-parity is defined 
by R p = (— i) 3B + L + 2S, ; S being the spin of a given particle, so that the particle content 
of the MSSM can be classified in SM particles with R p = +1 and superpartners with 
Rp = —1. The L- and B violating Yukawa interactions in the superpotential being 
odd under R p vanish. However, phenomenologically it is already sufficient to forbid 
the 5-violating interactions and to constrain the L-violating terms such that they are 

Note that the ep — > fi + jet + X events observed by HI and discussed in Jl^] fail significantly the 
kinematical constraints that should hold for the process eq — > LQ — > fiq. 
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Figure 10. Illustrative example 
of resonant squark production fol- 
lowed by i? p -violating and conserv- 
ing decays into e + and jets. 




consistent with the remaining, much weaker experimental bounds. 

In such a scenario one can have the term X' i:jk LiQjDk in the superpotential, where 
Li and Qj are the lepton and quark S , f/(2) i -doublet superfields, respectively, D k is 
the down-quark S77(2)£-singlet superfield, and A^ fc denote the dimensionless Yukawa 
couplings. The indices i,j, k refer to the fermion generations. Most interestingly 
at HERA, the _R p -violating couplings with i = 1 can give rise to resonant squark 
production p|J 0| In e + p scattering, the possible production channels are e + Uj — > d^ 
and e + dk — > uij- In particular valence d quarks allow to probe uij squarks and A^ 
couplings. The present indirect bounds on and Mg p6[ leave an interesting search 
window for HERA with the exception of A' m in which case the limit from neutrinoless (5(5 
decay precludes the observation of effects at HERA. The possible squark interpretations 
of the observed excess of events in NC DIS reported in ffl, H have been studied in detail 
in 



27, 28, 2S, 11 



In the following, we consider scenarios with a single nonvanishing coupling A^. 
Furthermore, we assume that gluinos are heavier than squarks, so that q — > g + q is 
kinematically forbidden. The resonantly produced squark can then decay into e + q 
via the ^,-coupling A^ resembling a scalar leptoquark, and into a neutralino Xi or 
a chargino xt pl us a quark via i? p -conserving gauge interactions. In the latter cases, 
the final state is determined by the subsequent neutralino and chargino decays. The 
heavier gauginos may cascade down to the lightest one, while the lightest neutralino 
may finally decay into e ± qq' or uqq, involving again the coupling A^. Similar |? p - decay 
channels exist for the x + with the exception of the e~-mode. An illustrative example 
is shown in Fig. [U]. In [ 23[ , a classification is given of all possible final states in terms 
of several distinguishable event topologies. Recently, the SUSYGEN ]30] package has 
been extended to include ^p-SUSY processes at HERA. The new version which will be 
released shortly allows a complete event generation for all possible channels. 



Further studies of ^p-processes at HERA can be found in (23[ |i[ |5) . 
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Figure 11. In- 
variant mass spectrum 
of e + multijet final 
states for data (symbols) 
and NC DIS expectation 
(histogram) . 
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4-2. Analysis 

The HI Collaboration |21[ analysed final states containing an and several jets. Since 
the lepton appears late in the cascade decay chain, it is expected to be degraded in 
energy. Consequently, the main selection criteria require an e ± candidate reconstructed 
at high y e together with two forward jets. The invariant mass spectrum of the 289 events 
selected is shown in Fig. 11 and compared with the distribution of the 285. 7± 28.0 events 
expected in the SM. The agreement is good. In the "wrong sign" channel, e + +p — > e _ + 
multijet, one candidate is observed, while 0.49 ± 0.2 are expected. 



4-3. Results 

Combining the LQ-like channel e + + jet with the e ± + multijet channels, one can put 
mass-dependent limits on the couplings X' 1 j 1 . To this end, four different sets of values 
have been chosen for the MSSM parameters M 2 and tan/3, leading to four typical 
scenarios for the lightest neutralino : photino-like states with mass 40, 100, and 150 GeV, 
and a zino-like state with mass 40 GeV. Given A' lj:l and the squark mass, the relevant 
widths and branching ratios are fixed. The HI analysis shown in Fig. [T2| excludes squark 
masses up to 262 GeV for A^ ~ 0.3. As expected, in the case of A' m , the HI bound 
is not competitive with the best indirect constraint coming from the non-observation of 
neutrinoless (3(3 decay. However, for A' 121 and A' 131 , the best indirect limits from atomic 
parity violation are improved by the HI direct bounds in a large region of the parameter 
space. Especially, for large Xi masses, the upper bounds on A' 121 and A' 131 are decreased 
by as much factor 4. 

The ftp squark scenarios considered above are also constrained by the first 
generation LQ searches at the TeVatron. However, these bounds become rather weak 
when (3{q — > e + q) is small, a possibility which exists due to the presence of R p - 
conserving decays. Recently, the TeVatron experiments have performed a dedicated 
Ijt p search |3l], [32] analysing final states with like-sign electrons plus jets. These may 
originate from the decay of pair-produced squarks. Assuming five degenerate squark 
flavours, squark masses below 252 GeV are excluded. If one squark (e.g. the lightest 
stop) is much lighter than the others, the cross section is significantly smaller. The 
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H1 Preliminary 



Figure 12. 95% 

CL upper limits on the 
coupling as a func- 
tion of the squark mass, 
for different masses and 
composition of the light- 
est neutralino Xv The 
most stringent indirect 
limits from neutrinoless 
(3(3 decay and atomic 
parity violation are also 
indicated. 
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corresponding mass bound is estimated [21J to decrease to about 150 GeV. 

5. Conclusions 

Results on neutral and charged current Deep Inelastic Scattering at high momentum 
transfer derived from the e + p data taken by the HI and ZEUS experiments in 1994- 
1997, have been presented and compared to the Standard Model expectations. The 
differential cross sections in terms of Q 2 are found to be in good agreement with the 
SM predictions. A slight excess at high Q 2 observed in the 1994-1996 data still persists 
with the 1997 data included, but with less significance. 

Using the same data, a search for leptoquarks and for squarks with R-parity 
violating couplings has been performed. No evidence was found for the existence of 
such species. The resulting bounds on leptoquark/squark masses and Yukawa couplings 
to lepton-quark pairs improve the constraints from other colliders and from low energy 
experiments in important regions of the parameter space. 
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